Formally exact quantization condition for nonrelativistic quantum systems 
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Based on the standard transfer matrix, a formally exact quantization condition for arbitrary 
potentials, which outflanks and unifies the historical approaches, is derived. It can be used to find 
the exact bound-state energy eigenvalues of the quantum system without solving an equation of 
motion for the system wave functions. 
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In general an exact solution for many quantum me- 
chanical problem is unavoidable and one is forced to 
resort to some types of approximation technique. The 
Wentzel-Kramers-Brillouin(WKB) approximation, based 
on the wave function expansion in powers of h, is stud- 
ied and applied frequently in many fields. Apparently, it 
only remains valid for the short-waves and then plays a 
role of a bridge between classical and quantum mechan- 
ics. Neglecting of the higher order of h directly leads to 
failing to describe the wave function at the classical turn- 
ing points, together with inaccurate energy eigenvalues 
in the lower bound states. LangerQ introduced an ap- 
proximation based on Bessel function, coming across the 
same problems as the WKB approximation. Rather than 
in terms of the exponential function, Miller'5| proposed 
a method by proper choice of the arbitrary equation. 
FriedrichQ believed the phase shift at the tuning point 
is energy-dependent with nonintegral Maslov indices in 
the limit of the long-waves. Recently, the structure of su- 
persymmetric quantum mechanics motivates a modified 
semiclassical quantization condition for one-dimensional 
Hamiltonians. But the supersymmetric WKB(SWKB) 
method only reproduces the exact energy eigenvalues for 
translationally shape-invariant potentials whose ground- 
state wave function is analytically known0 j^. When 
the two or three-dimensional potentials are dealt with, 
Langer-like corrections is needed, which is exact only 
for solvable spherically symmetric potentials. All above 
methods are unable to give an exact quantization condi- 
tion for arbitrary nonrelativistic quantum systems. 

In this paper, from a different starting point, we first 
take the standard transfer matrix as the basis of analy- 
sis, and deduce the exact quantization condition in which 
the scattering-led phase shift recovers the neglecting of 
series of higher order h in the WKB approximation. How- 
ever, the scattering-led phase shift term is of discrete sum 
formalism, furthermore, it is a bit purer mathematical 
and lacks of physical insights. But the defined equiva- 
lent attenuation vector is recognized to be nothing more 
than —ip (x) /ip (x) , so that an elegant integral formal- 
ism of the scattering-led phase shift can be instead of 
the discrete sum formalism. A Riccati equation is ob- 
tained, which is associated the wave function with the 
scattering-led phase shift, and we show the quantization 



condition outflanks and unifies the historical approaches, 
leading directly to the exact eigenvalues and the exact 
eigenfunction so that we can define an exact real quan- 
tum action-angle variable J (E) determining the energy 
eigenvalues. 

Without loss of generality, we assume a potential with 
two classical turning points xl and x^ which are deter- 
mined from V (xl) — V {xji) = E. Far from the two 
points, we truncate the potential at xc and x^. In or- 
der to employ the transfer matrix, [xc^^n] is divided 
into I + f + g equal section layers with width /i, letting 
[xc,xl] have I section layers, [xl, xr] have / section lay- 
ers and [xr^xd] have g section layers. In each section 
layer, the potential can be regarded as constant. Hence 
the transfer matrix of j th layer can be written as 



cos{Kjh) ^-^sm{Kjh) 
KjSm{Kjh) cos{Kjh) 



(1) 



where kj is defined as fiKj ~ ^ 2m \E — Vj\ and the sub- 
script J is an integer ranging from 1 to l + f+g. Vj stands 
for the potential energy at the j th section layer, m is the 
particle mass and E is the energy eigenvalue. Using the 
boundary condition that the wave function tjj (x) and its 
first derivative ip' (x) are continuous at the boundary be- 
tween two neighboring section layers, the matrix equation 
is given 



ip (xc) 



i+f+a 

n 



Ip [xd) 



(2) 



The right- and left-hand side in Eq.(2) are simultaneously 



multiplied by the matrix 



and divided 



by Ip (xc) V' {xd), then Eq.(2) is changed to the form 



i+f+g 

n 



= 0. 



(3) 



The following analysis is based on the section layer 
width /i — > 0, the mass m — I and Planck's con- 
stant h — 1. We define Pc = —tp'{xc)/4'{xc) 
and Pd = tp' {x d) / ip {x d) ■ Making use of the trans- 
fer characteristic of matrix in Eq.(3), the values of 
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—■ijj'(x)/4>{x) are calculated at any position x. After 
algebraic manipulations Eq.(3) is equivalent to 



w=l+l 

Pi 



= nvr + arctan 



arctan 



.1^1+1 

where 5 (n) is determined ft'om 

Pi+l 



l+f-1 

i=l + l 



arctan 



— arctan 



Kl+f 



Pi+l 



(4) 



with 



arctan | ''^^ ) ^ Kjh 



(5) 



(6) 



and Pj are subject to the boundary condition Pq = 
Pc, Pi+f+g+i — Pd- From Eq.(5), we know that if the 
potential is nonconstant, 6 (n) is nonzero, so that S (n) 
can be interpreted as the scattering-led phase shift which 
is in general energy-dependent or quantum number n- 



dependent. It is clear that ki+i 



0, while Pi and 



and Ki+f = ^2 [E-Vi+f+i] 
Pi+f+i is positive and finite for bound states, so 
arctan = arctan (P/_|_/_|_i/Pi+/) = 7r/2 that is 

the energy-independent phase shift at the turning points. 
If a wave is reflected at where V {x^) oo, mean- 
while. Pi +O0 and Ki+i is a finite positive number, so 
the half phase shift at Xl arctan {Pi/m+i) is also exactly 
equal to the same tt/2 as that at the turning point. The 
first term in Eq.(4) can be integrated as '^wh = 

J^^ ndx. Given that Eq.(6) is equivalent to Kjh — 
arctan (Pj+i/^j) — arctan (Pj/Kj), and the deriva- 
tive of arctan (P/k), namely, d [arctan (P/k)] /dx = 
{ndP - Pdn) / (k^ + p2V we have 



arctan 



arctan 



P, 



AP/t,-P,AK 

J , o2 ■ (7) 



Due to Ak 

reduces to a first order differential equation 

dP{x) 



J - = and AP = P^+i - P,-, Eq.(7) 
ferential eq 

^2 {x) + P2 {x) , (8) 



dx 



where 



P(x) 



tp{x) ' 



(9) 



which is a well-known Riccati equation. Substitu- 
tion Eq.(9) into Eq.(8) results in the time- independent 
Schrodinger equation ■0 (x) + k^i/j (x) = 0. From an- 
other angle, the deducing of Eq.(8) is proven to be true. 



Likewise, Eq.(5) can be integrated as 
, , dn P 

S{n) 



dx P2 + 

■ (10) 

where k — dn/dx and P = dP/dx. From Eq.(9), we 
obtain the wave function 



{x) = Nq exp 



P \ X ] dx 



(11) 



which is acquired to be square integrable. In order to 
determine P{x), one must solve Eq.(8). The initial con- 
dition of Eq.(8) is determined from Pd, while the func- 
tion ip (x) = Ni exp [—a (a: — xd)] is matched to the wave 
function decaying exponentially in the range [xd , oo] , 
with a — y^2m [V (xd) — E] which is the equivalent at- 

tenuated vector, so Pu = —ip (xd) /ip (xd) = ct- 

Now from Eq.(4) and Eq.(lO), i.e., J!^^ ndx + 5 (n) = 
{n + 1) TT, we have derived the exact quantization condi- 
tion 



k'-^ j da; = 2 (n-t- l)7r,n = 0, 1,2..., (12) 



where the integral is counterclockwise around a closed 
contour C which encloses the two turning points. The 
quantization condition complies with the requiring the 
total phase during one period of oscillation to be an 
integral multiple of 27r for a particle oscillating in the 
classically allowed region between the two turning points 
[a^L, xn]. Comparing with scattering theory in which the 
sub- waves phase is shifted when the particle is scattered 
by a potential, we can gain the better understanding why 
there is a term 5 [n) in Eq.(12), which is ignored by other 
approximations such as the WKB approximation and the 
Bohr- Sommerf eld quantization condition. 

In contrast with the classical Hamilton- Jacobi the- 
ory, the quantum action-angle variable eigenvalues J 
should be defined as J = J {E) = {1 / 2tt) <f^ K dx with 
K = K — K P/ P . Since — K P/P has the same mo- 
mentum dimension as k, K is referred to as the com- 
plete quantum momentum function (QMF). The point is 
emphasized that the J {E) is defined at arbitrary ener- 
gies, not only the eigenvalues. The wave function outside 
the classically allowed region, corresponding to whether 
the energy eigenvalues or other arbitrary energies, can- 
not be exponential growth, on the contrary, is bound to 
decay exponentially because the potential V{x) is larger 
than the particle E outside the classically allowed region, 
the particle cannot fully enter the classically forbidden 
region, implying the wave function should decay expo- 
nentially. As a result, there is no the points where log 
derivative P diverges. Considering the exponentially de- 
caying wave function outside the classically allowed re- 
gion and with the help of the Riccati Eq.(8), J{E) in 
Eq.(12) can be calculated for arbitrary E, clearly the be- 
havior of J{E) is monotonic in all E. The new definition 
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of J is a conceptual breakthrough not only because no 
one reached it from the founding of quantum mechanics 
to now but also Eq.(12) possesses the basis going back to 
the Bohr-Sommerfeld quantization condition in the old 
quantum theory. Thereby, a one-to-one correspondence 
between the invariant torus and quantum eigenvalues is 
unambiguously established. 

It should be noted that our defined J (E) is a real func- 
tion which describes the system motion, and, in addition, 
enables one to complete Hamilton's original program of 
associating a wave with a particle motion. In essence, 
the scattering of a particle in the bound state alters the 
momentum of the particle, as a result, contributes sig- 
nificantly to the exactness of Eq.(12). At the same time, 
Eq.(12) indicates the energy eigenvalue is involved both 
with the particle behavior between the two turning points 
and with the particle behavior outside the two turning 
points, which is the different aspect between the classi- 
cal and quantum mechanics, while the condition J (E) = 
n + 1 yields the exact energy eigenvalues, the number 1 
comes from the contribution due to the phase shift at the 
turning points. In view of the Einstein-Brillouin-Keller 
(EBK) quantization condition^ 01 since it only includes 
the particle behavior between the two turning points, it 
is impossible to generally give birth to the exact results. 

It is well worth mentioning the Bohm's hidden variable 
theory 0. A pair of equations is obtained, in which one of 
terms is called quantum potential, while we define a QMF 
K in Eq.(12). Apparently whether the so-called quan- 
tum potential or the QMF is related to the wavefunc- 
tion. According to Bohm's theory Newton's second law 

can be written as mcPv/dt^ — — V iV + O') \^ (*}[lot- 

where V is classical potential and Q is the quantum po- 
tential. But in our theory the classical k, and quantum 
momentum —k P/P , namely, the QMF contribute 
to the quantum dynamics of the system described by 
mdr/dt = k — n P/P \r=r{t)- This detailed investigation 
will appear elsewhere. 

The quantization condition Eq.(12) together with 
Eq.(8)-Eq.(ll) is exact both for the solvable potentials 
and for the unsolvable potentials, in order to illustrate 
them explicitly, let us first examine the detailed values of 
the scattering-led phase shift and detailed expressions of 
the phase integral for several familiar solvable potentials. 
Why these potentials are chosen is that their analytical 
wave functions are known so that one can calculate S (n) , 
consequently, one can confirm the exactness of Eq.(8)- 
Eq.(12). 

(1). One-dimensional infinite square well V (x) = at 
X C [0, L] but oo at other x. Due to dn/dx = 0, one has 
d (n) — and the quantization condition reads 

[ V2Edx ^ {n + l)Tr,n = 0,1,2..., (13) 
Jo 

which yields exact energy eigenvalues En = (n + 1)'^ Eq 



n 




-^present 


-^Maslov 





-0.97815416 


-0.97815416 


-0.97834291 


1 


-0.93556613 


-0.93556613 


-0.93566866 


2 


-0.87203511 


-0.87203511 


-0.87210568 


3 


-0.78795362 


-0.78795362 


-0.78800723 


4 


-0.68386490 


-0.68386490 


-0.68390852 


5 


-0.56051533 


-0.56051533 


-0.56055291 


6 


-0.41901295 


-0.41901295 


-0.41904703 


7 


-0.26131274 


-0.26131274 


-0.09251698 


8 


-0.09248716 


-0.09248716 


-0.09251698 



TABLE L Comparisons of exact, present, and the Maslov 
index energy eigenvalues for the radial Woods-Saxon potential 
^ ('') = - i+cxp [2f,-30)1 +^^ ' = 1' setting m = h = 1. 



with Eq = 7r^?i^/2L^. However, the WKB approximation 
\/2Edx = (n -I- 1/2) 7r?i can not do so. 
(2). One-dimensional harmonic oscillator V {x) = 
x^/2, whose eigenfunction can be expressed by ele- 
mentary functions, i.e., -00 = -^j'K^^'^ ex:^\—c?x^ 
■01 = \/2a/7r^/^aa; exp [— a^x^/2] substituting -0^ into 
Eq.(8) and Eq.(lO), we find 6 (n) = 7r/2, to a certain de- 
gree, which shows Eq.(lO) is exact because Eq.(12) re- 
duces to the familiar WKB quantization condition that 
remains exact for the one-dimensional harmonic oscil- 
lator. For the three-dimensional harmonic oscillator 
Vg^y (r) = r2/2 /2r2, substituting the known 

wave functions into Eq.(8) and Eq.(lO), we have 8 (n, I) = 
{21 1] 



2^T{iTi) 

tion reads 



/4, and the quantization condi- 



1)-{21-1) 



(n-t-l)7r, (14) 



V2m [E- Veffir)] in Eq.(14). 
the total Maslov index|3l /i happens to 



which yields the exact energy eigenvalues E = 2n + I + 
3/2. K is defined by k 
For this case, 

accord with the term 21 - 2y/r(r+ 1) + 3. (3). One- 
dimensional Coulomb potential V (x) = — l/|x|. Though 
there is a singularity at a; = 11], in the similar way, we 
find S {n) = TT, and the quantization condition reads 



ndx — riTT, n = 1, 2, 3. 



(15) 



where xr = —1/E. Eq.(14) gives the exact energy 
eigenvalues En = — 1/ (2n^). For the three-dimensional 
Coulomb potential Veff (r) = —l/r + l(l + l) /2r^, we 

find S {n, I) = 



condition reads 



-1/r + l{l + l)/2r^, 
TT, and the quantization 



ndr 



+ TT = {n+l)TT,n = 0,1,2..., 

(16) 



4 



E 



SWKB 



14.312470 



14.468822 
42.227424 
43.899629 
68.213256 
74.975448 
94.567829 

7 108.464359 

8 124.328691 



14.312470 
14.468822 
42.227424 
43.899629 
68.213256 
74.975448 
94.567829 
108.464359 
124.328691 



14.312670 
27.174940 
39.499957 
47.940588 
62.315504 
76.652144 
92.325304 
108.985251 
126.433812 



TABLE II: Comparisons of exact, present, and SWKB energy 
eigenvalues for the double oscillator V {x) = 10 (la;| — 3)^, set- 
ting m = Ti = 1. 



which determines the exact energy eigenvalues E = 
-1/ {2N^) with N = l + n+l. 

Eqs.(13)-(16) are novel and immediate expressions un- 
der no other auxiliary modifications, whose exactness 
can make one convinced that the scattering-led phase 
shift 5 (n, I) does exist objectively, and it plays a cru- 
cial role in determining the energy eigenvalues. Be- 
cause of it, when the WKB approximation is applied 
to the three-dimensional potentials, Langer-correction 
is needed. Interestingly, if the number I [l + \) is 
instead of {l + l/2f in Eq.(14) and Eq.(16), they 
are naturally changed to Langer-correction formalism 
^2m[E-V,ff {T)]dr = (n + 1/2) tt with K// (r) = 

V ir) + {1 + l/2f /2r2 . Moreover, when the Coulomb- 
like potentials is tackled, the problem of singularity is 
naturally avoided. With regard to the unsolvable poten- 
tials, generally 5 {n, I) is not an analytical expression of n 
and I. Next the energy eigenvalues for the unsolvable ra- 
dial Woods-Saxon potential^3| and the unsolvable dou- 
ble oscillator potential^^l with Eq.(12) are tabulated in 
Table I and II, respectively, and the complete agreement 
with exact results can be seen, implying Eq.(12) is exact. 
-gexact obtained through numerical technique. Dealing 
with the double oscillator potential, as the energy of the 
particle is relatively smaller, four turning points can oc- 
cur, assuming xi < X2 < xy, < x^. Because the particle 
will be motioning in the range [a;i,a;4], the quantization 



condition J^^" (^n — n P / dx — {n + 1) -k gives the en- 
ergy eigenvalues. For other complex potentials with more 
then four turning points, in the similar way, the contour 
is taken from the smallest turning point to the biggest 
turning point. In these cases, the tunneling effect takes 
place. 

In summery, based on the transfer matrix rather than 
the wave function expansion in powers of Ti, we have de- 
rived the exact quantization condition for arbitrary non- 
relativistic quantum systems and have shown it unifies 
and transcends the historical approaches. Naturally, we 
have obtained the corresponding real quantum action- 
angle variable, contrasted with the classical Hamilton- 
Jacobi theory, which discloses the quantum-classical cor- 
respondence. 
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